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ABSTRACT: Linkage disequilibrium (LD) plays an important role in genomic selection and mapping quantitative trait loci (QTL). In 
this study, the pattern of LD and effective population size (/VJ were investigated in Chinese beef Simmental cattle. A total of 640 bulls 

were gcnotypcd with lUuminaBovinSNPSOBeadChip and lUuminaBovinHDBeadChip. We estimated LD for each autosomal 
chromosome at the distance between two random SNPs of <0 to 25 kb, 25 to 50 kb, 50 to 100 kb, 100 to 500 kb, 0.5 to 1 Mb, 1 to 5 Mb 
and 5 to 10 Mb. The mean values of were 0.30, 0.16 and 0.08, when the separation between SNPs ranged from 0 to 25 kb to 50 to 100 
kb and then to 0.5 to 1 Mb, respectively. The LD estimates decreased as the distance increased in SNP pairs, and increased with the 
increase of minor allelic frequency (MAF) and with the decrease of sample sizes. Estimates of effective population size for Chinese beef 
Simmental cattle decreased in the past generations and was 73 at five generations ago. (Key Words: Chinese Beef Simmental Cattle, 
Bovine Genome. Linkage Disequilibrium. Single Nucleotide Polymorphisms. Effective Population Size) 



INTRODUCTION 

Linkage disequilibrium (LD) denotes non-random 

association between alleles at different loci. LD is the 
theoretical basis of genomic selection (GS) and genome- 
wide association study (GWAS), that is also important in 
gene mapping, estimates for effective population size, 
population structure and so on (Nachman, 2002). Molecular 
markers such as single nucleotide polymorphisms (SNPs) 
and microsatellites were widely used to estimate the extent 
of LD. The level of LD is usually influenced by non-genetic 
factors and genetic factors containing genetic linkage, 
selection, the rate of recombination, the rate of mutation, 
genetic drift, non-random mating and population structure. 

The effective population size (A^^) is defined as the 
number of individuals in an ideal population that would 
show the same amount of dispersion of allele frequencies 
under random genetic drift or the amount of inbreeding as 
in the population under consideration, and is usually less 
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than the absolute population size (Wright, 1938). A^^is an 
important parameter, as it can help to explain how cattle 
populations evolve and expand, and by definition describe 
the rate of inbreeding accumulation and loss of genetic 
variation. Estimates for A^,, can be obtained from 
heterozygote excess or LD. Presently, estimates for A'^ based 
LD data are are more frequently used than heterozygote 
excess, and therefore complement evolutionary studies of 
cattle populations (Hayes et al., 2003). 

Recently the discovery of large numbers of SNP 
through sequencing of the cattle genome has generated 
extensive research in quantifying LD characteristics (Farnir 
et al., 2000; Odani et al., 2006; McKay et al., 2007; 
Sargolzaei et al., 2008; Kim and Kirkpatrick, 2009; Qanbari 
et al., 2010). A recent report showed that high density 
markers were used to study the extent of LD in Angus, 
Charolais and Crossbred beef cattle (Lu et al., 2012). 
However, similar studies were not reported in Simmental 
cattle which are an important economic breed of beef cattle. 
LD indicates population characteristics and has different 
pattern on each chromosome. Hence, it is necessary to study 
the extent of LD, and then to estimate effective population 
size in Simmental cattle. 

China's role in international beef markets has grown 
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significantly in the past years, and domestic production is 
projected to continue to increase (Longworth et al., 2001). 
However, China does not have special-purpose beef cattle. 
To increase beef production, American, Canadian and 
Australian Simmental cattle have been introduced into 
China and crossbred with native dual-purpose Simmental 
cattle, which are named Chinese beef Simmental cattle. In 
current research, high density SNPs data from Chinese beef 
Simmental cattle were used to analyze the pattern of LD, 
and to infer the effective population size up to 2000 
generations ago. Meanwhile, we evaluated the effects of 
minor alleUc frequency (MAF) and sample size on LD 
estimations. 

MATERIALS AND METHODS 

Animals 

Experimental animals consisted of 640 young 
Simmental bulls, bom in 2008 to 2010, originated from 
Ulgai, located at Xilingol league, Inner Mongolia, China. 
DNA was extracted from blood of the bulls using the 
routine procedures. The lUuminaBovineHD chip was used 
to genotype 504 young bulls and their autosomal 
chromosomes contained total of 735,293 SNPs. 
Additionally, 136 young buUs were genotyped with 
IlluminaBovineSNP50, and 51,582 SNPs were detected on 
their autosomal chromosomes. There were 46,000 SNPs in 
common between two chips. In the present study, quality 
control standards for SNPs data were Hardy-Weinberg 
equilibrium (p>10"^), MAF>0.05, SNP call rate >0.95 and 
Mendel error rate <0.05. 35,079 common SNPs survived 
after being filtered on quaUty control standards, which were 
used to analyze the extent of LD. 



and to present or process this information in various ways. 
In this study, we used the command plink -cow -bfile 
filename -Id-window -r2 0 -out outname. To display the 
decay of LD, distances of pair-wise SNPs were binned into 
seven types of intervals (0 to 25 kb, 25 to 50 kb, 50 to 100 
kb, 100 to 500 kb, 0.5 to 1 Mb, 1 to 5 Mb and 5 to 10 Mb) 
along the first 10 Mb of each chromosome, and mean / was 
computed for each interval. Table 2 shows information for 
all the SNP pair groups. 

Three factors, chromosomes, MAF and sample sizes, 
affecting LD estimation were studied based on / data 
computed above. 

Genetic distance 

In the high-density SNP chip, genetic distance for SNP 
pairs could not be obtained. Therefore, physical distance 
was used to replace genetic distance for the estimation of 
effective population size in the ciurent study. 100 kb of 
physical distance in genetic distance is approximate 
equivalent to 0. 1 cM. SNP physical position from the UMD 
3.1 bovine assembly (http://www.ncbi.nlm.nih.gov/ 
assembly/313678/) was used in this study. 

Effective population size estimation 

LD data make it feasible to estimate Ng. Sved (1971) has 
proposed the relationship formula for LD and as follows: 

_ 1 ^ 1 (Equation 2) 

kN^c + a n 

^12 (Equation 3) 

^2 In, kc 
(r — )kc 

n 



LD estimation 

Several statistics parameters were proposed to measure 
the extent of LD. D' (Lewontin, 1964) and (Hill, 1974) 
were widely used in practice, but their functions are 
different, was considered to be a better descriptor of LD 
as it is more robust and not sensitive to changing gene 
frequency and effective population size (Terwilliger et al., 
2002; Zhao et al., 2007). 

Assume two loci A and B, each locus has two alleles 
(denoted Ai, A2 and Bi, B2, respectively). Pai, Pa2> Psiand 
Pflzare the frequency of each of the alleles. Pn, P12, P2iand 
P22show the frequency of haplotypes AjBi, A1B2, A2Biand 
A2B2. Thus, / can be expressed as: 

^2 ^ (PiiPi2-Pi2P2i)^ (Equation 1) 



Where Ne is the effective population size, r is the mean 
value of LD estimation value of SNP pairs. A: is 4 or 2, 
which denotes autosomal chromosomes or sex 
chromosomes respectively, c represents the genetic distance 
in Morgan (M). Generation is calculated with l/2c. n is the 
chromosome experimental sample size. «= 1 (Sved, 1971) 
in the absence of mutation, otherwise a = 2 (Lewontin, 
1964; Weir and Hill, 1980; McVean, 2002). In the majority 
of studies, the formula for the absence of mutation is chosen 
to estimate Ne- Hence, in our study, k = A and c = 1 were 
chosen. Seven types of SNP pairs with the physical 
distances 25 kb, 50 kb, 100 kb, 500 kb, 1 Mb, 5 Mb and 10 
Mb were respectively chosen to estimate the A^^ of Chinese 
beef Simmental cattle from 5 generations ago. 

RESULTS 



PLINK (PurceU et al., 2007) includes a set of options to SNP statistics 
calculate pair-wise hnkage disequiUbrium between SNPs, SNPs information for every autosomal chromosome is 
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Table 1. Statistical information for analyzed SNP 


Chromosome 


Length 

(Mb) 


Number of SNP 


Average SNP Interval 

(Mb) 


Longest interval 

(Mb) 


Shortest interval 
(kb) 


BTAl 


158.14 


2,290 


0.05 


1.13 


0.13 


BTA2 


136.52 


1,815 


0.06 


1.47 


0.08 


BTA3 


121.37 


1,726 


0.05 


1.32 


0.11 


BTA4 


120.61 


1,731 


0.05 


0.59 


4.90 


BTA5 


119.73 


1,454 


0.06 


0.61 


0.15 


BTA6 


119.21 


1,769 


0.05 


1.60 


1.36 


BTA7 


112.61 


1,543 


0.05 


1.49 


0.00 


BTA8 


113.36 


1,629 


0.05 


0.52 


1.80 


BTA9 


105.46 


1,381 


0.06 


0.93 


0.45 


BTAIO 


104.21 


1,491 


0.05 


2.41 


0.28 


BTAll 


107.04 


1,514 


0.05 


1.10 


0.88 


BTA12 


91.09 


1,128 


0.06 


3.34 


0.24 


BTAl 3 


84.15 


1,242 


0.05 


1.95 


0.38 


BTAl 4 


84.61 


1,205 


0.05 


3.62 


0.11 


BTA15 


85.05 


1,141 


0.05 


0.86 


2.90 


BTAl 6 


80.92 


1,071 


0.05 


2.39 


0.17 


BTA17 


74.96 


1,080 


0.05 


0.81 


4.80 


BTAl 8 


65.97 


928 


0.05 


0.97 


5.60 


BTA19 


64.01 


967 


0.05 


0.54 


1.40 


BTA20 


71.79 


1,096 


0.05 


0.51 


3.50 


BTA21 


70.61 


921 


0.05 


1.22 


0.90 


t> 1/\ZZ 


oyj.D^ 


oOj 


yj.yjD 


u.^o 


u.zu 


BTA23 


52.13 


749 


0.05 


1.14 


4.20 


BTA24 


62.64 


910 


0.05 


0.47 


0.10 


BTA25 


42.80 


687 


0.04 


0.49 


1.30 


BTA26 


51.68 


732 


0.05 


0.73 


6.60 


BTA27 


45.36 


671 


0.05 


1.28 


11.00 


BTA28 


46.19 


657 


0.05 


2.14 


0.02 


BTA29 


50.97 


688 


0.05 


0.91 


1.80 



given in Table 1 . The total autosomal chromosome length of 
Chinese beef Simmental cattle was 2,54L30 Mb. The 
longest Bos taurus autosomal chromosome is BTAl (length 
= 158.14 Mb), and the shortest is BTA25 (length = 42.80 
Mb). 35,079 common SNPs between two chips covered the 
whole genome in this study. Average adjacent SNPs spacing 
was 54.17±61.44 kb, and the largest spacing situated on 
BTA14 was 3620 kb (between ARS-BFGL-NGS-37733 and 
Hapmap42739-BTA-95927). The mean MAF of the genome 



was 0.28+0.13, and followed an almost uniform distribution, 
as can be seen in Figure L 

Extent of LD across the genome 

The mean values of for each autosomal chromosome 
for distance bins of 0 to 25 kb, 25 to 50 kb, 50 to 100 kb, 
100 to 500 kb, 0.5 to 1 Mb, 1 to 5 Mb and 5 to 10 Mb were 
calculated. Table 2 shows that the average r^is 0.30, 0.23, 
0.16, 0.08, 0.05, 0.04 and 0.03 at different distance bins for 
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Figure 1. Minor allele frequency (MAF) distribution for total SNPs. 
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Table 2. Statistical information for LD over various distances 



Distance 



Number of SNP pairs 



Average 



0-25 kb 


4,100 


0.30 


25-50 kb 


12,412 


0.23 


50-100 kb 


18,855 


0.16 


100-500 kb 


92,378 


0.08 


0.5-1 Mb 


179,984 


0.05 


1-5 Mb 


368,467 


0.04 


5-10 Mb 


654,755 


0.03 



r : denotes the extent of LD. 

Simmental cattle, respectively. Figure 2 shows the LD 
decay over varying distances of the genome. The measured 
LD was high for pairs of SNPs within close proximity. 
However, there is a strong LD in the long distance SNP 
pairs. 

The extent of LD was significantly different among 
chromosomes. The average for SNPs separated by 
intervals 0 to 25 kb, 25 to 50 kb, 50 to 100 kb, 100 to 500 
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Figure 2. LD decay for Chinese beef Simmental cattle in whole 
autosomal chromosome. 

kb, 0.5 to 1 Mb, 1 to 5 Mb and 5 to 10 Mb in each 
autosomal chromosome are presented in Table 3. The mean 
value of for distances less than 25 kb was 0.30, but higher 
for BTA9 and BTA21 (0.363 and 0.364, respectively), and 
lower for BTA27 (0.209). The average of was 0.30 in 



Table 3. Statistical information for average r as distance between pairs of SNP up to 10Mb for the genome 



CHR 










SNP pairs Distance 








0-25 Kb 


25-50 Kb 


50-100 Kb 


100-500 Kb 


0.5-1 Mb 


1-5 Mb 


5-10 Mb 


BTAl 


0.32210.298 


0.25510.257 


0 


18610.201 


0.08310.103 


0.04910.034 


0.03910.022 


0.03410.015 


BTA2 


0.33510.317 


0.25110.252 


0 


16510.179 


0.08510.097 


0.05110.037 


0.04010.023 


0.03410.015 


BTA3 


0.30210.291 


0.22910.235 


0 


16610.188 


0.07810.089 


0.04810.032 


0.03810.023 


0.03310.016 


BTA4 


0.27710.267 


0.24210.245 


0 


17310.186 


0.08910.108 


0.05310.042 


0.04010.023 


0.03310.015 


BTA5 


0.28210.282 


0.25010.255 


0 


16910.197 


0.08510.109 


0.05510.045 


0.04410.029 


0.03510.018 


BTA6 


0.32210.318 


0.24910.257 


0 


17410.195 


0.08910.103 


0.05810.047 


0.04410.029 


0.03510.017 


BTA7 


0.31810.296 


0.25310.251 


0 


17910.189 


0.08510.106 


0.05110.041 


0.04110.029 


0.03410.016 


BTA8 


0.29110.287 


0.24610.256 


0 


15810.177 


0.07910.096 


0.04910.032 


0.03910.022 


0.03410.015 


BTA9 


0.36310.321 


0.23410.250 


0 


16610.188 


0.07310.083 


0.04810.032 


0.03910.023 


0.03410.015 


BTAIO 


0.28610.306 


0.23110.241 


0 


16410.187 


0.07210.082 


0.04510.031 


0.03710.018 


0.03310.014 


BTAU 


0.31810.291 


0.25110.257 


0 


16910.195 


0.07910.092 


0.05110.040 


0.04110.026 


0.03410.016 


BTA12 


0.25710.258 


0.20710.231 


0 


16010.188 


0.07410.086 


0.05010.036 


0.04010.024 


0.03410.016 


BTA13 


0.29110.306 


0.18410.196 


0 


13710.160 


0.06910.080 


0.04710.034 


0.03710.020 


0.03210.014 


BTAU 


0.29010.299 


0.23310.230 


0 


16210.181 


0.08110.099 


0.04910.035 


0.03810.021 


0.03310.014 


BTA15 


0.28510.293 


0.21310.229 


0 


16410.184 


0.07610.085 


0.04910.037 


0.03910.023 


0.03310.017 


BTAl 6 


0.33010.304 


0.24810.238 


0 


16510.190 


0.07510.087 


0.04810.034 


0.03710.019 


0.03310.014 


BTA17 


0.28910.288 


0.22710.241 


0 


15110.167 


0.07510.085 


0.04810.033 


0.03810.021 


0.03310.014 


BTA18 


0.34310.318 


0.23110.237 


0 


14010.169 


0.07310.077 


0.04810.030 


0.03910.021 


0.03310.013 


BTA19 


0.26010.245 


0.20110.222 


0 


14810.183 


0.06910.079 


0.04610.030 


0.03710.020 


0.03210.017 


BTA20 


0.23410.241 


0.21510.226 


0 


15210.175 


0.07410.081 


0.04910.034 


0.03910.023 


0.03410.018 


BTA21 


0.36410.308 


0.24810.231 


0 


17810.192 


0.07910.090 


0.05010.036 


0.04010.025 


0.03310.022 


BTA22 


0.33710.324 


0.22810.227 


0 


14310.169 


0.07410.081 


0.0510.036 


0.0410.022 


0.03310.014 


BTA23 


0.27710.272 


0.19610.200 


0 


14210.167 


0.07710.104 


0.05210.046 


0.03810.019 


0.03210.013 


BTA24 


0.32910.286 


0.25910.251 


0 


16310.191 


0.07610.080 


0.04910.031 


0.0410.023 


0.03410.015 


BTA25 


0.26510.250 


0.1910.220 


0 


12910.148 


0.06410.068 


0.04510.026 


0.03610.018 


0.03110.012 


BTA26 


0.29110.279 


0.23110.239 


0 


15410.183 


0.07010.071 


0.04910.035 


0.03910.021 


0.03310.017 


BTA27 


0.20910.220 


0.19310.221 


0 


13210.163 


0.06410.070 


0.04810.033 


0.03710.019 


0.032+0.013 


BTA28 


0.22510.253 


0.19810.225 


0 


1310.145 


0.06110.058 


0.04410.028 


0.03610.019 


0.03210.012 


BTA29 


0.29610.278 


0.20610.233 


0 


13010.152 


0.06410.072 


0.04610.030 


0.03810.023 


0.03210.013 



CHR denotes chromosome, r: Means+SE. 
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Figure 3. Average / estimates at different physical distances for three different minor allelic frequency (MAP) thresholds. Mean LD 
estimates are pooled over all autosomal chromosomes, and three different minimum threshold cut off levels for minimum allele 
frequency are shown. 



SNP pairs with physical distances of <25 kb and decreased 
to 0.16 at distances of 50 to 100 kb, this result was similar 
to that previously reported (Qanbari et al., 2010; Lu et al., 
2012). A similar study found the extent of LD (r = 0.59) in 
approximately 50 kb on north American Holstein cattle, 
which was much larger than that found in our study 
(Sargolzaei et al., 2008). 

MAF and LD 

In this study, three different minimum allelic frequency 
(MAF) thresholds (0.05, 0.1 and 0.2) were used to study the 
effects of MAF on the extent of LD. Figure 3 shows that 
MAF has a significantly effect on the mean value of r^, 
especially over short distances (0 to 25 kb). The mean 
values of increase significantly with an increasing MAF. 
For example, from 0 to 25 kb, the mean value of for 
MAF>0.05 was 0.24, however, with MAF>0.1 and 0.2, the 



mean value of increased to 0.29 and 0.34, respectively. 

Sample size and LD estimates 

As can be seen in Figure 4, sample sizes affect the LD 
estimation value. In this paper, five different sample sizes of 
n = 25, n = 50, n = 100, n = 200 and n = 400 were randomly 
selected from the total set to study the effect of sample size 
on estimates of the level of LD. The mean / were greater 
when sample size is smaller, and this phenomenon is more 
noticeable for LD estimation across a SNP interval more 
than 500 kb. There were no significant differences for LD 
estimates when sample sizes were greater than 400 and SNP 
distances less than 50 kb. 

Effective population size 

The extent of LD for different chromosome fragment 
length could reflect the effective population size of different 




0-25kb 25-50kb 50-lOOkb 100-500kb 0.5-lMb l-5Mb 
SNP pairs Distance (Mb) 



5-lOMb 



Figure 4. Average estimates at different physical distances for six different sample sizes. Mean LD estimates are pooled over all 
chromosomes, and six different sample sizes are shown. 
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Table 4. Statistical information for effective population sizes of Simmental cattle 

SNP pairs distance 

25 kb 50 kb 100 kb 500 kb IMb 5 Mb IMb 



Genetic distance 0.00025 0.0005 0.001 

Generations ago 2,000 1,000 500 

N, 2377 1697 1344 
Genetic distance: Morgan. JV« = Effective population size. 

past generations. Table 4 shows of Simmental cattle in 
past generations. Estimates of A^^, for 2,000 generations ago 
was approximately 2,377 and down to 73 at 5 generations 
ago. Estimates of for Chinese beef Simmental cattle 
show an increasing trend when plotted against increasing 
past generations (Figure 5). 

DISCUSSION 

Recent developments in high-throughput SNP panels 
have generated enthusiasm and interesting in GS and 
GWAS on cattle. Linkage disequilibrium maps can increase 
power and precision in association mapping. Qanbari et al. 
(2010) reported an average level LD of 0.30 over pair wise 
distances less than 25 kb based on 40,854 SNPs in 810 
German Holstein cattle. Kim and Kirkpatrick (2009) 
reported LD of >0.80 over genomic regions of 
approximately 50 kb using 7119 SNPs in North America 
Holstein cattle. Lu et al. (2012) reported the extent of LD in 
Angus, Charolais and crossbred beef cattle based on 
lUumina BovineSNP50_v2 Beadchip and lUumina 
BovineSNP50_vl Beadchip, with the level of LD being 
0.29, 0.22 and 0.15 when the distance range between 
markers is 0 to 30 kb, respectively. This could be attributed 
in part to the difference in populations between the current 
study and previously reported research. Furthermore, in the 
current study, we used 35,079 SNPs distributed across the 
entire bovine autosomal chromosome for the analysis of LD 



0.005 0.01 0.05 0.1 

100 50 10 5 

611 484 123 73 



in Chinese beef Simmental cattle. The statistic denotes 
the extent of LD. The extent of LD showed a decreasing 
tendency with increasing distances of the genome. The 
mean was much higher between close loci, and the result 
was the same as previously reported estimates (Farnir et al., 
2000; Smith et al., 2006; Kim and Kirkpatrick, 2009; 
Qanbari et al., 2010; Lu et al., 2012). However, a low level 
of LD can exist between two SNPs that are closely adjacent, 
while markers that are more distant can show a higher than 
expected level of LD. This situation also appeared in 
Unkage disequilibrium studies on human and model animals 
(Reich et al., 2001). It could be caused by selection, the rate 
of recombination, mutation and genetic drift (Nachman, 
2002). 

The mean values were different for the same fragment 
length on different autosomal chromosomes. Higher LD 
was found for BTA21. This may reflect selection for traits 
that are strongly influenced by QTL on this chromosome in 
this breed. Chinese beef Simmental cattle are a popular 
breed in Chinese beef production and genetic trends suggest 
a strong selection for growth and meat traits. A majority of 
studies have shown highly significant evidence for the 
presence of QTLs affecting meat traits (McClure et al., 
2010) on BTA21. In addition, when selection operates at a 
locus, the neighboring loci in close linkage with the locus 
under selection will have an enhanced extend of LD. When 
selection occurs at multiple loci in epistasis, LD between 
loci under epistatic selection and their tightly linked loci 
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Figure 5. Estimated Ng for Chinese beef Simmental cattle over time from linkage disequilibrium data. 
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will be created and enhanced (Du et al., 2007). 

Estimates of LD across the whole genome could be 
affected by many factors. In this study, removing SNPs with 
very low MAFs also lead to lower numbers of SNPs 
available for study, which can also lead to bias of LD 
estimates. There are several published papers observing a 
similar phenomenon in other species (Khatkar et al., 2008; 
Yan et al., 2009; Qanbari et al., 2010). LD estimation is 
greater with MAF increasing at a short SNP pairs distance, 
but the phenomenon is not sensitive when SNP pairs 
achieve a distance of 1 Mb. Sample size is another factor 
that affects estimation of the extent of LD (Khatkar et al., 
2008; Yan et al., 2009). A small sample size (n = 25) can 
also lead to the biased estimates for LD. However, there are 
no significant differences for the mean when sample sizes 
exceed 100, especially when the given extent interval of LD 
is less than 100 kb. In addition, previous research on 
Holstein cattle demonstrated that a sample of 400 or more 
was required for reUable estimation of LD (Khatkar et al., 
2008). A similar study in humans that found sample sizes 
would be even higher, which may be due to humans having 
a larger effective population size (Chen et al., 2006). 

Hill (1974) proposed a method for estimating effective 
population sizes. In this method, estimates of A^,, depend on 
the number of animals alive at any time and the variance of 
progeny number per sire. In addition, previous research 
showed that the latter played a key role in the decrease of 
the population size (Mukai et al., 1989; Nomura et al., 
2001). To miiximize the net response in economic merit for 
dairy cattle, FAO (1998) reported an effective population 
size of 50 per generation was required to maintain the 
fitness in a breed. Goddard and Smith (1990) suggested a 
minimum effective number of 10 bull sires per generation, 
equivalent to 40 individuals per generation. McParland et al. 
(2007) used this traditional method to estimate Ne for 
550,591 Ireland Simmental cattle, the result showed that 
was 127 at the current generation. However, pedigree 
information was often missing or error that caused the 
decUne of A^^ estimated accuracy. In our study, the estimate 
for N^, of Chinese beef Simmental cattle was approximately 
73 for 5 generations ago, well above the reported numbers. 
This could be attributed to a sufficiently large number of 
sires being used to produce animals in the current dataset, 
and thus a small variance of family size was generated. The 
slope of the A^^ suggests that the population sizes were 
decreasing consistently fast, possibly due to the use of 
artificial selection, and therefore actions is required to 
maintain a larger Ng. 
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